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Abstract. ESO 4m class telescope and VLT deep imaging of the isolated neutron star RX J0720. 4-3125 reveals 
a proper motion of /i = 97 ± 12mas/yr and a blue U-B color index. We show that a neutron star atmosphere 
model modified to account for a limited amount of hydrogen on the star's surface can well represent both the 
optical and X-ray data without invoking any additional components. The large proper motion almost completely 
excludes the possibility that accretion from interstellar medium is the powering mechanism of the X-ray emission. 
It also implies that the proposed spin down is entirely due to magnetic dipole losses. RX J0720. 4-3125 is thus a 
very likely middle aged cooling neutron star. Its overall properties are quite similar to some of the long period 
radio pulsars recently discovered, giving further support to the idea that RX J0720. 4-3125 may be a pulsar whose 
narrow radio beam does not cross the Earth. 



1. Introduction 

ROSAT has identified a small group of hot and radio-quiet 
isolated neutron stars (INSs, see reviews by Motch [2001 
Treves et al. EUUTl and Haberl l2UU3|) . None of these objects 
seems associated with a SNR and low interstellar absorp- 
tion towards them indicates distances of the order of a few 
hundred parsecs at most. Grating X-ray spectra of the two 
brightest members, RX J0720. 4-3125 (Paerels et al-HED 
and RXJ1856. 5-3754 (Burwitz et al. I2UUT1 Drake et al. 
2002}, do not significantly depart from pure blackbody en- 
ergy distributions. Inferred blackbody temperatures range 
between 50 and 100 eV. Soft X-ray pulsations with periods 
around 10 s are now detected in a majority of sources al- 
though some period candidates still require confirmation. 
Optical counterparts are very faint objects whose bright- 
ness seems to be mostly due to thermal emission from the 
neutron star surface (see Pavlov et al. 120021 for a recent 
review). These INSs undergo little, if any, magnetospheric 
activity and are thus prime targets to measure fundamen- 
tal parameters through modeling the neutron star thermal 
emission and to constrain the debated equations of state. 

The evolutionary status of these objects is still unclear. 
Some of the 10 s to 10 9 INSs born in the Galaxy could 
have space velocities and rotation periods slow enough 
to accrete from relatively dense phases of the interstellar 
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medium and be re-heated to temperatures similar to those 
observed (see e.g. Blaes & Madau IT993fl . The fraction of 
old INSs detectable in X-ray is of the order of several 
percents but strongly depends on the assumed velocity 
distribution (Madau & Blaes 1994). This hypothesis cur- 
rently faces two main difficulties. First, the space density 
of hot radio-quiet INSs is much lower than predicted by 
the accretion models (Treves et al. l2001fl . Second, the ro- 
tation periods, although relatively long for neutron star 
standards, are still too short to allow accretion to take 
place in average interstellar medium conditions unless the 
magnetic fields of INSs have significantly decreased by 
about two orders of magnitudes over their life time (see 
e.g. Wang lLW7|) . 

Alternatively, these INSs could be young cooling ob- 
jects. Long rotation periods are observed in few radio pul- 
sars as well as in AXPs showing that under certain cir- 
cumstances (efficient braking or birth with a low angu- 
lar momentum) , relatively young neutron stars can rotate 
slowly. 

Optical observations can efficiently help to discrimi- 
nate between the two hypotheses, re-heating by accretion 
from the insterstellar medium or relatively young cool- 
ing objects. For some of the objects, proper motions can 
be measured which, assuming Bondi-Hoyle accretion, may 
rule out the accreting scenario. Additionally, the optical 
to X-ray flux ratio is sensitive to the chemical composition 
of the neutron star surface (Pavlov et al. I1996J) and may 
therefore reveal the presence of newly accreted material 
from the interstellar medium or from a fossil disc. 
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In both the optical and X-ray band, RX J0720. 4-3125 
is the second brightest INS among those discovered with 
ROSAT. This object exhibits pulsations at a period of 
8.39s (Haberl et al. ITWI) with |P| < 3.6 lCT^ss" 1 
(Kaplan et al. I2002[l . The X-ray spectrum is blackbody- 
like with kT°° ~ 86 eV (as given by the ROSAT data). 
ESO-NTT and Keck observations have led to the optical 
identification of RX J0720. 4-3125 with a faint stellar-like 
object (Motch & Haberl 1998. Kulkarni & van Kerkwijk 
1998|.We report here on new 4m-class telescope and ESO- 
VLT observations of RX J0720. 4-3125 which reveal its 
proper motion, show its blue UV continuum and provide 
constraints on the nature of its thermal emission. While 
this paper was in the refereeing process we became aware 
of a preprint by Kaplan et al. H2003[l reporting on HST 
observations of RXJ0720. 4-3125 which confirm the steep 
optical to UV energy distribution of the star) . 

2. Optical observations 

Optical imaging was obtained over five ESO periods span- 
ning the time interval from February 1997 till January 
2003. The log of observations is given in Tabled Data col- 
lected in 1997 have been already presented and discussed 
in Motch & Haberl ifTMHjl . In January 1999 we used the 
ESO-NTT and the SUSI2 imager equipped with EEV44- 
80 CCDs #45 and #46 (2048x4096 pixels). The January 
2000 run was carried out at the ESO 3.6m using EFOSC2 
and the LORAL 2048 2 CCD #40. At UT1 we used in both 
cases FORS1 and a TK 2048 2 CCD. Individual exposures 
of - 10 mn were obtained in 1997, 1999, 2000 and 2002- 
2003 for the B and R bands while the U band 2000 images 
were 15 mn long. In January 1999 and 2000, the U and B 
filters were alternated and in all instances the telescope 
was moved by a few arcseconds in different directions be- 
tween two exposures. Raw images were corrected for bias 
using over-scan regions or mean bias frames and corrected 
for flat-field using dithered sky images obtained at dawn. 
Individual corrected images were then moved to a common 
frame using a set of several nonsaturated bright stars and 
then stacked together applying a statistical criterion in or- 
der to reject cosmic-ray impacts. Because of the excellent 
seeing prevailing during ESO-VLT observations, the geo- 
metrical distortion effects of the FORS camera (up to 0.2" 
on the edge of the field) combined with dithering slightly 
deteriorated the quality of the stacked images. Distortion 
effects were thus corrected before stacking using the midas 
rectify image task and parameters taken from the FORS 
users Manual (Issue 1.5). After correction, the positional 
errors of reference stars dropped from 30 mas to 5 mas and 
the FWHM image size sharpened by 0.1 to 0.05". 

3. Optical photometry 

Fig. ^ shows the deep B image obtained by summing 
the 3.4 h FORS1 observation in 2002-2003 of best image 
quality resulting in an average FWHM seeing of 0.56". 
The optical counterpart of RX J0720. 4-3125 (object "XI" 



from Motch & Haberl ifTMSjl or "X" in Kulkarni & van 
Kerkwijk (1998)) is clearly detected, thanks to the excel- 
lent images and VLT throughput. Our relatively deep U 
exposure obtained by summing the 7.2 h long observation 
in January 2000 collected at the ESO-3.6m (see Fig. 0) 
confirms the strong UV excess of the counterpart seen by 
Kaplan et al. lj3rM|> . Object "X2" in Motch & Haberl 
(199$J) does not appear exactly at the same position in 
the 1997 image as in the Keck observation of Kulkarni & 
van Kerkwijk (1998) or as in our 1999, 2000 January and 
VLT images. We believe that a low amplitude unfiltered 
sky background fluctuation could be the origin of its dis- 
placement. The revised position of "X2" outside of the 
ROSAT error box and its lower UV flux excludes it as a 
possible counterpart of RX J0720. 4-3125 (Kulkarni & van 
Kerkwiik ll998|l . Our R image in January 1999 is not deep 
enough to detect the counterpart. 

Absolute photometric calibration was obtained in the 
course of the January 2000 run. Standard stars in the fields 
of PG 0231+051, PG0942-029 and PG1047+003 (Landolt 
1992) were used to derive U and B magnitudes of three 
comparison stars named A, B and D located close to 
RXJ0720. 4-3125. Stars A and B are the same as in Haberl 
et al. I|1997ll . Their magnitudes are listed in Tableland 
their positions are shown in Fig. ^ The sky background 
close to RX J0720. 4-3125 is somewhat contaminated by 
stray light from a nearby group of relatively bright stars. 
Because of its equatorial mounting images obtained at 
the ESO 3.6m further suffer from fixed diffraction spikes 
caused by the secondary spider. Following Kulkarni & van 
Kerkwijk l|1998|l . we found that the most accurate pho- 
tometric measurements were achieved by fitting a 2-D 
Gaussian distribution of fixed width on top of an inclined 
plane representing the sky background and stray light 
from nearby stars. The fixed width of the Gaussian profile 
was derived from bright nonsaturated field stars and the 
photometric scale measured on averaged images since the 
flux from bright stars is not accurately conserved in the 
image stacking process. Colour transformations provided 
by ESO were systematically applied in order to compen- 
sate for the large colour differences between RXJ0720.4- 
3125 and the photometric field reference stars. Magnitudes 
derived from the various runs are listed in Table [3] The 
1999-2003 B magnitude of RX J0720. 4-3125 is consistent 
with that reported by Kulkarni & van Kerkwijk ( 1998), B 
= 26.6 ± 0.2, and the re-analysis of the 1997 images using 
2-D Gaussian fitting instead of aperture photometry also 
provides consistent values, B = 26.58 ± 0.25. The U-B 
index is very blue, U-B = -1.11 ± 0.18, assuming no B 
optical flux variation. The simultaneous but less accurate 
B measurement in January 2000 yields a consistent U-B 
= -0.76 ± 0.23. Our estimate of the U band flux agrees 
well with the optical to far UV distribution measured by 
Kaplan et al. (|2003|l . The difference of 0.17 magnitude in 
B between the two ESO-VLT observations is only a 2.6 a 
effect and is thus probably not significant. We also used 
the high quality optical data from the last ESO-VLT run 
to search for night to night variability. There again, we 
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Table 1. Optical data 
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Table 2. Magnitudes of field reference stars (la errors) 

Star U B 

A 20.585 ± 0.014 20.542 ± 0.009 
B 21.011 ± 0.015 20.841 ± 0.009 
D 18.770 ± 0.014 18.766 ± 0.008 



Table 3. U and B magnitudes of RX J0720.4-3125 (la 

errors) 

Epoch U B 

February- April 1997 26.58 ± 0.25 

January 1999 26.79 ± 0.20 

January 2000 25.68 ± 0.17 26.44 ± 0.15 

December 2000 26.787 ± 0.040 

Dec 2002 - Jan 2003 26.620 ± 0.050 



find some scatter with a maximum amplitude of 0.25 mag 
from one night to the other which is only significant to 
a 2.9 a level. It is more likely that the rather large pho- 
tometric scatter is due to a slight underestimate of the 
actual photometric errors. 

4. Proper motion study 

We searched for possible proper motion by comparing the 
positions of the optical counterpart at five epochs from 
1997 till 2003. Apart from the 1997 images which are the 
sum of data collected over a time interval of 2 months 
from February 7 till April 4, the four other epoch obser- 
vations took place during integer number of consecutive 
years within 27 days from the same calendar date. Overall, 
the change in position due to trigonometric parallax is ex- 
pected to be of the order of 16 mas at maximum for an as- 
sumed distance of 100 pc (mainly due to the spread in time 
of the 1997 observations), and can therefore be neglected 
for our purpose. As for the photometry we found that 
fitting 2-D Gaussian profiles was providing the most ac- 




Fig. 1. Sum of best seeing B band images obtained with 
FORS1 on UT1 in December 2002 and January 2003. 
Positions of the photometric reference stars are shown. 
North is to the top and East is to the left. The arrow 
marks the position of RXJ0720. 4-3125 



curate results and good estimates of the positional errors. 
Positions from automated processes such as Sextractor or 
DAOPHOT turned out to be less reliable in general. The 
fixed width of the stellar profiles in right ascension and 
declination were first derived by averaging the values ob- 
tained on a set of relatively bright but nonsaturated stars. 
Then the position of a group of 10 reference stars was mea- 
sured, together with that of RX J0720. 4-3125 and of 17 
faint test objects. Reference stars are evenly spread within 
30" of RXJ0720. 4-3125 and are in the range of B magni- 
tude from 21.8 to 18.8. Comparison objects are distributed 
in the same sky region as the reference stars and were se- 
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Table 4. Maximum errors (milli-arcseconds) in RA and 
Dec for the astrometric reference stars 
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Fig. 2. The summed U band image obtained with 
EFOSC2 on the ESO-3.6m in January 2000. The area cov- 
ered is roughly the same as that seen in Fig. 1. The arrow 
shows the position of RX J0720. 4-3125 

lected for their stellar-like appearance. Their magnitudes 
range from B ~ 27.0 to 23.5. Using a rotation, translation 
and scaling transformation, the 1997, 1999, 2000 January 
and 2003 positions were then moved to the 2000 December 
FORS1/VLT reference frame which provides the most ac- 
curate positions. Not all the reference stars are common to 
all frames and the 1997 to December 2000 transformation 
was computed with only 8 stars. We checked that the qual- 
ity of the transformation was not depending heavily on the 
number of astrometric reference stars used. Since part of 
the proper motion could be due to bulk motion of the ref- 
erence stars we also computed the astrometric solution of 
the December 2000 VLT image using 78 USNO-A2 field 
stars. Reference stars do not show any systematic differ- 
ence in position in the time interval from 1980 (USNO-A2) 
till 2000. The accuracy with which reference stars coincide 
after transformation is rather good. Their maximum offset 
to mean positions after transformation are listed in Table 
H Between December 2000 and 2002-2003 VLT observa- 
tions, the mean residual of the 10 reference stars is only 
4.5 mas. We plot in Fig. [31 the apparent displacements 
of the 17 faint comparison stars together with those of 
RX J0720. 4-3125 and of the astrometric references derived 
from the two FORS1/VLT observations. None of the ref- 
erence stars exhibit significant motion with respect to the 
rest of the group. Compared to stars of similar brightness, 
the apparent motion of RXJ0720. 4-3125 is clearly of high 
significance. 

We also corrected for differential refraction ef- 
fects resulting from the large colour difference between 
RXJ0720.4-3125 (U-B = -1.11) and the reference astro- 
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Fig. 3. Apparent proper motion of local astrometric refer- 
ence stars (open squares), comparison objects (stars) and 
of RXJ0720.4-3125 (filled square) using December 2000 
and December 2002 - January 2003 VLT images. 



metric stars (mean U-B = 0.17). In the B filter, this trans- 
lates into a shift of ~ 130A in effective wavelength. Since 
all observations were done at low airmass, the corrections 
remain small, of the order of 5 to 30 mas and are usu- 
ally negligible compared to positioning errors. Differential 
refraction effects are only significant in the case of the 
December 2000 VLT observation. 

We then fitted a proper motion solution expressed in 
terms of total motion and direction angle to the vari- 
ous positions of RX J0720. 4-3125 using a minimum x 2 
algorithm. The best fit shown in Fig. 01 is obtained for 
fj, = 97 ± 12 mas/yr. The direction of motion is oriented 
29° ± 7° north of west (all errors are 90% confidence level 
for one parameter). 

5. Optical to X-ray energy distribution 

In the X-ray domain RXJ0720. 4-3125 was observed with 
the ROSAT, SAX, XMM-Newton and Chandra observa- 
tories. As results inferred from different instruments are 
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Fig. 4. Positions of RXJ0720. 4-3125 for the five epochs 
(filled squares) together with the best fit proper motion 
solution (open squares) corresponding to /i = 97 mas/yr 

subject to calibration uncertainties (see, e. g., Burwitz et 
al. EED for the case of RXJ1856. 5-3754), we used only 
the Chandra data collected with the LETGS grating in- 
strument in a total exposure of 37.5 ks in February 2000. 
First results from the spectral analysis of these data have 
been presented in Pavlov et al. I|2002|l . They found that 
neither light-element (hydrogen or helium) nor heavy- 
element (pure iron or solar mixture) neutron star atmo- 
sphere models can explain the multi-wavelength spectral 
distribution of the emission detected from RXJ0720.4- 
3125. The net result of this analysis was remarkably sim- 
ilar to the situation with the spectral fitting of the op- 
tical and X-ray data on RX J1856. 5-3754: while provid- 
ing acceptable fits to the X-ray data the light-element at- 
mosphere models yield unrealistic large radii of the INS 
and significantly over-predict the actual optical fluxes (by 
a factor of ~ 300), whereas the heavy-element models 
do not fit the X-ray data due to numerous features in 
the model spectra caused by absorption at atomic transi- 
tions in ions of the heavy elements (see Zavlin & Pavlov 
12002 for details). Pavlov et al. I^0Cj2) found that the best 
fit to the Chandra LETGS data on RX J0720. 4-3125 is 
provided by a blackbody model with the temperature 
kTg^ = 80 — 82 eV emitted from an area of a radius 
RJ^ = (2.0 - 2.2) (D/lOOpc) km (D is the distance to 
the source), absorbed with the hydrogen column density 
Nh = (1.5 - 1.7) x 10 20 cm" 2 . 

Similar to the case of RX J1856. 5-3754, the best black- 
body model derived from the X-ray data on RX J0720.4- 
3125 underpredicts the optical fluxes detected from the 
source by a factor of 4. Pavlov et al. (2002) specu- 
lated that a two-blackbody model, originally suggested 



by Pons et al. (2002) to describe the multi-wavelength 
data on RXJ1856. 5-3754, may be also applied to the ra- 
diation from RX J0720. 4-3125, with the "soft" component 

of kT bb, s < 43 eV and R bb, s > 6- 1 (D/100pc) km. Such 
a model, with X-rays originating from a hot area on the 
star's surface and the optical fluxes emitted from the rest 
of the surface, may explain the pulsations of the X-ray ra- 
diation from RX J0720. 4-3125. In addition to the thermal 
components, HST observations by Kaplan et al. (2003) 
reveal evidence for a nonthermal power law component. 

However, while this simple multiple-component model 
seems to be in agreement with the properties of the emis- 
sion detected from RX J0720. 4-3125, it is hardly reconciled 
with the fact that, as stars are not blackbodies, radiation 
emitted by a star should deviate from a blackbody model. 
Also, in the case of RX J1856. 5-3754, the absence of pulsa- 
tions put severe constraints on the two-blackbody model, 
requiring either a particular geometrical configuration or 
strong gravitational deflection (Ransom et al. 2002 ) . 

All previous models of neutron star surface radiation 
were based on a conventional assumption that there is 
enough matter (e. g., pure hydrogen or iron, or a mixture 
of elements) on the neutron star surface to make the at- 
mosphere layers optically thick at all energies of interest. 
A typical estimate for such an amount in terms of the 
total surface column density is ytot > 10 — 100 g cm -2 
(depending on the surface temperature) to provide the 
equilibrium (or diffusion) solution of the radiative trans- 
fer problem in the very deep layers (see Mihalas I1978fl . 
Under this assumption, spectra of the emitted radiation 
are solely determined by the temperature distribution in 
the atmosphere (which grows towards larger depths) and 
do not depend on properties of star's layers lying under- 
neath the atmosphere. In the case of a light-element at- 
mosphere composition, the model spectra were found to 
be much harder at higher photon energies than blackbody 
ones providing the same radiative flux. The reason for this 
effect is that high-energy photons with longer mean-free- 
paths are emitted from deep surface layers with tempera- 
tures larger than the so-called effective temperature, the 
fourth power of which determines the total energy flux 
(see, e.g., Zavlin et al. and Zavlin & Pavlov E003 

for details). Due to this property the model spectra, when 
applied to observational data, usually yield lower temper- 
atures by a factor of 2 — 3 and much larger emitting areas 
by a factor of 50 — 200 than estimates obtained from black- 
body fits. However, it is not known a priori how much, for 
example, hydrogen has been deposited on the surface of 
a neutron star (e. g., by accretion from the interstellar 
medium or from a fallback disk). If the hydrogen layer is 
not optically thick at all energies, but only at lower ones, 
then the atmosphere structure (temperature run) is mod- 
ified, that in turn should affect the spectra of emergent 
radiation. 

We investigated nonmagnetic atmosphere models with 
various values of the total column density of hydrogen y to t 
in the atmosphere. We assumed the atmosphere to be in 
the radiative and hydrostatic equilibrium, and that the 
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radiation at the inner boundary y=ytot is given by the 
diffusion solution (although, the latter condition is not 
strictly justified — see Discussion). The main result we 
obtain is rather natural: reducing the total column den- 
sity lowers the temperature of surface layers which emit 
high-energy photons, that in turn leads to the softening of 
the Wien tail in the model spectra, making them close to 
the blackbody ones with the same effective temperature. 
This means that optical fluxes predicted by these models 
applied to X-ray data on INSs are expected to be signif- 
icantly decreased in comparison with those given by the 
standard model, but still larger than values yielded by the 
blackbody model. We computed hydrogen models in a grid 
of effective temperature and total column density (down to 
v tot = 10~ 3 g cm" 2 ) and applied these models to the UV 
- optical and Chandra LETGS data on RX J0720. 4-3125 
taking into account the nonthermal power law discovered 
by Kaplan et al. (20Q2J). A best fitting model is shown in 
Fig. El The parameters of the model are: total atmosphere 
column density yt t = 0.16 g cm -2 , effective temperature 
kT ff = 57 eV, distance to the star D = 204 pc (assuming 
the standard neutron star radius R = 10 km and mass 
M = 1.4M Q ), and interstellar absorption N H = 1.1 X 10 20 
cm~ 2 . As seen in Fig. [SJ the model nicely fits the LETGS 
data (x 2 =0.94 in the X-ray domain) and agrees well with 
the detected optical fluxes. The total mass of hydrogen 
obtained in the fit is 4 7rR 2 y tot = 2xl0 12 g. Assuming a 
stellar age of 10 6 yr, this translates into a mean accretion 
rate of ~ 0.6 g/s, which is many orders of magnitude lower 
than predicted by the accretion model. 

6. Discussion 

Despite the fact that the thin atmosphere model provides 
a good fit to the observational data in the broad energy 
range at reasonable neutron star parameters, it serves 
rather as an illustrative example because of a number of 
shortcomings in the simplified approach it has been de- 
rived with. First of all, the assumption on the equilibrium 
solution at the maximal depth in such an atmosphere may 
not be justified at higher X-ray energies where the emit- 
ted radiation could be affected by properties of the stellar 
surface under the atmospheric layer. It means that the 
boundary condition at the inner boundary y= y to t , taken 
currently to be the same as in the standard modeling (with 
much larger values of ytot), may need to be modified. For 
instance, one could consider such a thin atmosphere on top 
of a solid neutron star surface or a thicker layer composed 
of another chemical element. Second, the model should be 
developed with account for the presence of a strong surface 
magnetic field, B ~ 10 13 G, as recently suggested by Zane 
et al. 1)2002(1 . It is safe to assume that in the case of mag- 
netized models, reducing the atmosphere thickness would 
also soften the spectrum of the emergent radiation. Since 
in strong magnetic fields the emergent spectra of the stan- 
dard models are intrinsically softer than the non-magnetic 
ones (Zavlin & Pavlov 2002), one may expect that the 
magnetized model needs a larger atmosphere thickness in 
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Fig. 5. Multiwavelength energy distribution of 
RXJ0720.4-3125. The optical and far UV data are 
from Kulkarni & van Kerkwijk (1998), Kaplan ct al. 
(2003) and this paper. The X-ray spectrum is that from 
the Chandra LEGTS data. Solid line shows the best thin 
neutron star atmosphere (NSA) + nonthermal power law 
(Kaplan et al. 2003, dashed line) model fitting both the 
X-ray and optical data (see text for the parameters) . The 
dotted line gives unabsorbed fluxes. 



order to describe the broad-band properties of the emis- 
sion from RXJ0720. 4-3125. 

Another point here is that this model is constructed 
assuming a uniform neutron star surface, which makes 
the model emission unpulsed, whereas the observed X- 
rays from RX J0720. 4-3125 reveal pulsations with a pulsed 
fraction of about 10%. A way to model the pulsations is 
to consider a nonuniform temperature distribution over 
the surface (see, e.g., Cropper et al. I2D0TI since such a 
temperature nonuniformity can be induced by the strong 
magnetic field (Greenstein & Hartke H983(l . Finally, in its 
present state, the thin atmosphere model applied to the 
multiwavelength data on RXJ1856. 5-3754 underpredicts 
the distance to this object by a factor of 2. In light of 
all this, the presented spectrum should be considered as 
a first, perhaps rather crude step towards more advanced 
modeling of radiation from a neutron star surface covered 
with a thin atmosphere layer, which might help to eluci- 
date the origin of thermal emission detected from INSs. 
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A proper motion of 97 ± 12 mas per year corresponds 
to a transverse velocity of Vt ~ 50 x (D/lOOpc) kms -1 . 
A small part of this displacement could be due to solar mo- 
tion (about 12 mas/yr in the right ascension if the source 
is at a distance of 100 pc). RX J0720. 4-3125 is thus the 
second radio-quiet soft X-ray bright isolated neutron star 
for which proper motion is detected. The value reported 
here is about a third of that measured for RX J1856. 5-3754 
(Neubauser Um Walter I200T1 Kaplan et al. EEE2J . 

Such a spatial motion, although modest for neutron 
star standards, most probably rules out the accretion sce- 
nario for RX J0720. 4-3125. The observed bolometric lu- 
minosity of 2.6 x 10 31 (D/100pc) 2 erg s" 1 (Haberl et al. 
I1997fl implies accretion rates of the order of M ~ 1.2 10 11 
(D/100pc) 2 gs _1 . On the other hand, the accretion rate 
can be estimated assuming the Bondi-Hoyle model: M ~ 
n v^Q 10 11 gs^ 1 where n is the mean density of the in- 
terstellar medium (in cm -3 ) and vio the relative velocity 
in units of 10 kms" 1 (Treves et al. I2001fl . Therefore, very 
high ISM densities of the order of n ~ 150 (D/lOOpc) 
cm -3 are required to account for the observed X-ray flux 
in the accretion scenario. This density is much higher than 
the mean density derived from soft X-ray photoelectric 
absorption or from optical reddening of field stars, n ~ 
0.2 - 0.4 cm" 3 (Motch & Haberl ITOOgjl . However, there 
have been recent reports on evidence for dense ISM struc- 
tures (n ~ 10 - 1000 cm -3 ) on small scales ranging from 
10 to 10 6 AU in regions of otherwise low mean densities 
(Lauroesch & Mever 120001 and Meyer & Lauroesch lT09"0|) . 
RXJ0720. 4-3125 could be presently passing through one 
of these high density structures at a velocity of > 10 x 
(D/lOOpc) AU yr -1 . Since such high density clouds are 
likely to be rare, we would only detect X-ray emission from 
old neutron stars when they pass through these small re- 
gions. The X-ray ionization of these cloud-lets may be de- 
tectable in the optical. It is not clear, however, whether the 
apparent absence of X-ray flux and temperature changes 
from RXJ0720. 4-3125 is compatible with the size of these 
dense structures and the expected time scale for accretion 
instability. This scenario furthermore, cannot account for 
the large P > 10 _14 ss _1 measured by Zane et al. ( 2002). 

Another implication of the high proper motion is 
that in no case can the spin-down derived by Zane et 
al. I|2()()2(l be due to a propeller torque extracted from 
the interstellar medium. Following Colpi et al. ©98), 
propeller spin-down can be written as P ~ 3.1 x 

10 -16 n 9/13 v "27/13 £8/13 p2 i /13 g g _i where ^ j g ^ pQ _ 

lar magnetic field in units of 10 12 G. At P = 5.4 10" 14 s s" 1 
the propeller mechanism dominates over the magnetic 
dipole losses if vio < 0.9 n 1 / 3 . Since vio is likely greater 
than unity, the propeller mechanism could only dominate 
if the star presently goes through a high density medium 
with n > 100. 

The possibility that RX J0720. 4-3125 breaks through 
propeller interaction with a fallback disk seems as well 
unlikely. First, assuming that the entire X-ray luminosity 



is due to friction leads to much larger P than observed 
(Zane et al. 120021 Kaplan et al. 12002) . 

In the framework of the blackbody interpretation of 
the X-ray data, a part of the optical flux excess above the 
extrapolation of this model to the low-energy band (see 
§ 5 and Fig. [5) could arise from the reprocessing of X- 
rays in a residual debris disk when the neutron star has 
entered the pulsar phase. Emission spectra of such discs 
were computed by Perna et al. H2000JI in a number of con- 
figurations. Although in the B band, the entire flux excess 
above the X-ray blackbody extrapolation could be due to 
a remnant disc, its energy distribution is rather red and 
would fail to account for the U band excess. Therefore, it 
is not clear whether such discs can significantly contribute 
to the overall optical emission of RX J0720. 4-3125. Similar 
conclusions are reached by Kaplan et al. ( 2003 ) . 

The lack of parallax and radial velocity as well as the 
large uncertainty on the parameters of the proper motion 
do not allow to guess the birth place of RX J0720. 4-3125 
as it is apparently possible for RX J1856. 5-3754 (Walter & 
Lattimer 2002|. The star is moving away from the Galactic 
plane but is still very close to it. Its general direction of 
motion is consistent with a birth either in the Sco OB2 
complex, as for RX J1856. 5-3754 (more probably in the 
lower Centaurus Crux region at D ~ 170 pc), or in the 
Vela OB2 + Trumpler 10 association (D= 360 - 410 pc). 
In both cases, the flight time from birth place is of the 
order of 10 6 yr with a rather large uncertainty. This value 
is generally consistent with the spin down age and with 
cooling time (see Zane et al. 120021 . 

Several high magnetic field (B > 10 13 G) radio pulsars 
with long periods have been recently discovered (Camilo 
et al. 120001 Morris et al. 2002). One of these radio pul- 
sars, PSR J1830-1135 has properties astonishingly close to 
those of RX J0720. 4-3125 (P = 6.22 s, P = 4.7xl0" 14 ss- 1 
and B = 1.7 x 10 13 G, Morris et al. 12002)1 . It could 
thus well be that RX J0720. 4-3125 and several other X-ray 
bright and radio-quiet INSs are radio pulsars whose radio 
beams do not cross the Earth (Motch 120011 Kaplan et al. 
12001 Kaplan et al. 120031 ). The effect of radio beam nar- 
rowing with increasing periods (Biggs 1990) could plausi- 
bly explain the much smaller ratio of long-period to short- 
period radio pulsars compared to that found for X-ray 
pulsars. 

7. Conclusions 

ESO-VLT and 4m class telescope imaging reveal that 
RXJ0720. 4-3125 has a proper motion of fi = 97 ± 12 
mas/yr. Such a relatively high spatial velocity almost com- 
pletely rules out the possibility that the neutron star ac- 
cretes matter from the interstellar medium. This strongly 
favours the idea that RXJ0720. 4-3125 is a cooling INS of 
an age around 10 6 yr. A propeller mechanism on a fallback 
disk, or on accreted material from the ISM cannot account 
for the braking detected by Zane et al. 1)2002(1 and thus 
the observed spin-down is likely due to magnetic dipole 
losses. Recent discoveries of radio pulsars with compara- 
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ble periods and magnetic fields combined with the effect of 
the radio beam narrowing at longer periods suggest that 
RX J0720. 4-3125 could well be an off-beam radio pulsar. 
As in RX J1856. 5-3754, the optical continuum lies well 
above the extrapolation of the blackbody model seen in 
X-rays. To account for the optical emission one can add 
a lower temperature blackbody component cold enough 
not to be detected in the X-ray range (Pavlov et al. 2002 
Kaplan et al. 12003(1 . However, neutron stars are unlikely 
to radiate like blackbodies. Spectra predicted by heavy- 
element atmosphere models have a number of prominent 
features which should be conspicuous in the Chandra and 
XMM-Newton data obtained with high-energy resolution. 
On the other hand, the standard (thick) hydrogen atmo- 
sphere models overpredict the optical emission by several 
magnitudes. We propose here a new kind of model atmo- 
sphere in which the thickness of the light-element layer has 
reduced values. Although the model still has a number of 
shortcomings, it can adequately fit the overall energy dis- 
tribution with a single component, requiring the distance 
to RXJ0720.4-3125 to be about 200 pc which still has 
to be compared with the actual value yielded by future 
parallax measurements. 
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